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ABSTRACT 



We have obtained a deep i-band image of the blazar 3C 279 while the target was in a low optical state. Due to the faintness of the 
optical nucleus we have made the first detection of the host galaxy. The host galaxy has an apparent I-band magnitude of 18.4±0.3 
and an efi'ective radius of (2.7 ±1.1) arcsec. The luminosity of the host galaxy Mf, = -23.8 is consistent with the luminosities of 
other radio-loud quasar host galaxies. Using the empirical correlation between bulge luminosity and central black hole mass Mw, we 
estimate log(M/v,/Me) = 8.9 ± 0.5, broadly consistent with values obtained by photoionization methods. 
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1. Introduction 

The host galaxies of active galactic nuclei (AGN) provide im- 
portant information of the environment and evolution of their 
central black holes. The close connection between the cen- 
tral black holes and their host galaxies is clearly manifested 
by the correlations between the central black hole mass 
and the host galaxy parameters, such as the central velocity 
dispersion cr ( Ferrarese & MerrittI I2OO0I: iGebh ardt et ar'2000*: 
Merritt & Ferrarese 2001), bulge mass Mbuige (Magorrian et al. 
19981) and galaxy light concentration Crf (1/3) dGraham et alj 
200 lb . It is usually much easier to determine the host galaxy 



parameters th an to apply kine matic methods or reverberation 
mapping (e.g. lKaspi et al.l l2000') making the host galaxy method 
an appealing alternative to estimate black hole masses in galax- 
ies. Host galaxies can also be used to study the unified schemes 
of AGN since their prop erties do not depend on viewing angle 
( IUrrv&Padovanilll995l) . 

3C 279 (z = 0.538) is an AGN which belongs to the class 
■ of blazars, whose most notable features are strong continuum 
' radiation and high variability over the whole electromagnetic 
spectrum and high optical and radio polarization. The blazar 
class comprises of flat spectrum radio quasars (FSRQs) and BL 
Lacertae objects (BL Lacs) whose main difference in the optical 
lies in the strength of their broad emission lines: the latter group 
has very weak or nonexistent emission lines whereas in the for- 
mer they are more pronounced. 3C 279 is classified as a FSRQ, 
although both the emission line strengths and y-ray properties 
place it close to the F SRQ-BL Lac border dKoraflcar et al.lll998t 
iGhiselUni et 10120091) . 

Despite of 3C 279 being subject to intensive study over the 
last three decades, determining its host galaxy parameters has 
received little attention. The only pu blished host galaxy study 
was made by iKotilainen et alj (Il998l) who observed 3C 279 in 
the H-band, but could only derive an upper limit of H > 13.7 
for the host galaxy. One of the main factors affecting the host 
galaxy detection efficiency is the nucleus to host ratio: when the 
nucleus is very bright compared to the host galaxy its is difficult 



to detect the host galaxy, let alone reliably determine its main 
parameters (magnitude and effective radius). 3C 279 is known to 
be highly variable in the optical b ands with historica l variations 
over 6 magnitudes in the B-band jWebb et al.l ll990). The rapid 
variability of 3C 279 makes the host galaxy detection even more 
challenging as is it difficult to observe the target in a low state. 

We are currently monitoring over 40 blazars in the R-band in 
the Tuoria Observatory blazar monitoring progranfl- In the end 
of June 2008 we observed 3C 279 to go into a fairly deep optical 
minimum (R ~ 16.8) and initiated prompt i-band imaging at the 
Nordic Optical Telescope (NOT) to detect its host galaxy. The 
results of this imaging are presented in this research note. 

Throughout this paper we use the cosmology //q - 70 km 
s ' Mpc ', Om = 0.3 and Qa = 0.7. 



2. Observations and data reduction 

The observations were made in June 23-24, 2008 with the Nordic 
Optical Telescope (NOT) using the ALFOSC instrument. This 
instrument has a pixel scale of (y.'189/pixel and a gain of 0.726 
e /ADU and a readout noise of 3.2 electrons. The field of view 
of the instrument is 6.5' x 6.5'. Altogether 24 images of the 
3C 279 field were obtained through an i-band filter with almost 
an uniform transmission between 725 and 875 nm. Individual 
exposure times were kept low enough (100-200s) to leave 3C 
279 and a nearby star (star 1 in Fig.[Tli unsaturated. The position 
of 3C 279 on the CCD was changed between exposures to be 
able to make a fringe correction image. 

The individual images were first bias-subtracted and then 
flat-fielded with twilight flats in the usual way using IRAlfl The 
fringe pattern visible in the images was coiTected by subtract- 
ing a fringe correction image produced by median combining all 



http ://users . utu.fi/~kani/ 1 m/index. html 
- IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Table 1. The results of the model fitting. 





Model 




Ihost 


'cff 






Nucleus 


15.99 






1.38 




Nucleus + De V. 


16.13 


18.43 


2'.' 7 


1.03 




Nucleus + disk 


16.08 


19.04 


2'.'8 


1.05 



Fig. 1. The summed i-band image of the field around 3C 279. 
The field size is 5.3'x5.4', North is up and East is to the left. 
Stars 1-4 are discussed in the text. 



24 images and applying simultaneously a sigma clipping cut to 
the pixel intensity values. After correction for the fringe pattern, 
the individual frames were registered and summed. The summed 
image shown in Fig. [T]has a FWHM of 0.79 arcsec and a total 
exposure time of 3300 s. 

Calibrati on of the fiel d was obtained via star 1 in Fig.IT] for 
which iSmith & Balonekl (199^) give I = 15.00 + 0.04. Since 
the i-band filter used in our observations closely matches the 
Cousins I-band filter, we expect only small color effects between 
our i-band magnitudes and standard I-band magnitudes. We con- 
servatively estimate the error of calibration to be 0. 1 mag. 



3. Analysis and results 

We looked for the host galaxy of 3C 279 by fitting two- 
dimensional surface brightness models to the observed i-band 
image. We used three different models: an unresolved nucleus 
only, nucleus -i- a de Vaucouleurs profile host galaxy (Sersic = 
0.25) and nucleus + a disk host g alaxy (P - 1.0). Det ails of this 
fitting process can be found in Nilsson et al.l (Il999h . In short, 
the model has three adjustable parameters, the magnitude of the 
nucleus /core, the magnitude of the host galaxy /host and the ef- 
fective radius of the host galaxy rgff . The unresolved nucleus is 
assumed to be centered on the host galaxy. The three parameters 
are adjusted via an iterative Levenberg-Marquardt loop until the 
minimum chi squared between the model and the observed sur- 
face brightness distribution is found. The fit was extended to an 
outer radius of 9'.'5, except that the pixels affected by a galaxy 
5'.'9 SE of 3C 279 were masked out from the fit. 

Prior to computing the chi squared the model was convolved 
with the PSF. We used star 1 in Fig. [1] for the PSF du e to its 
high si gnal to noise and proximity to 3C 279. As in Nilsso n et alj 
( l2008l) . we studied the variability of the PSF across the FOV 
by extracting the surface brightness profiles of stars 1-4 in Fig. 



[T] and deriving the rms scatter between the profiles. The PSF 
variability was then parameterized by a parabolic expression 



crpsp(r) = 0.035 + 2x 10" 



(1) 



where crpsF('") is the uncertainty in the PSF relative to the inten- 
sity at radius r (pixels) from the PSF center. We then computed 
the expected variance cr- in a pixel with intensity / (ADU) and 
distance r from the center of 3C 279 from the expression 



cr2 = 



(2) 



where G is the effective gain and R is the effective readout noise, 
and compute the from 



X' 



(3) 



where M is the model intensity and the summation is over all 
unmasked pixels within the fitting radius. By including the PSF 
variability into the computation of we can ensure that the 
results are not dominated by PSF errors, which are most pro- 
nounced close to the center of the object. 

The results of this model fitting are summarized in Table [T] 
and Fig. |2] We see a clear excess in 3C 279 over the PSF and 
this excess is clearly above any PSF variability across the FOV 
(see the lower panel of Fig. |2]i. The models with a host galaxy 
give a much better fit than a model with pure nucleus with x^ 
indicating a near-perfect fit. The model with a de Vaucouleurs 
host galaxy gives a slightly better fit that the model with a disk 
galaxy but the difference is not significant according to our er- 
ror simulations (see below). However, since the de Vaucouleurs 
profile gave formally a better fit and no blazar host galaxy has 
ever been associated with a disk galaxy, we concentrate in the 
following to the results with the de Vaucouleurs profile. 

We have studied the sensitivity of the results to random 
noise, PSF variability and errors in the assumed profile of the 
host galaxy by performing 50 model fits to simulated images of 
3C 279. In these simulations /core, 4o,st and r^s were held con- 
stant at the values shown in Table [T] but the contribution from 
the different noise sources changed from one simulation to an- 
other 

The random noise in these simulations is assumed to raise 
from readout and photon noise according to the first term on 
the right hand side of Eq. |2] To model the PSF variability we 
created two PSF models in each simulation, slightly differing 
from each other The simulated model was convolved with the 
first PSF and the model fits were made with the second PSF. Both 
PSFs consisted of an elliptical Moffat profile with (3 - 2.5. The 
ellipticity of the first PSF was randomly drawn from an uniform 
distribution with a minimum of 0.0 and a maximum of 0. 12 and 
the position angle from an uniform distribution between and 
180 degrees. For the second PSF the ellipticity was drawn from 
a Gaussian distribution with a mean equal to the ellipticity of the 
first PSF and cr - 0.05. The position angle of the second PSF 
was similarly drawn from a Gaussian distribution with a mean 
equal to the first PSF and cr = 12.5 degrees. The limits and 
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^ , , , Table 2. Previous black hole mass determinations for 3C 279. 




I I I I I I ">-- I 

0123456789 




-0.2 I 1 1 1 1 1 1 1 1 I 

0123456789 

radius (arcsec) 



Fig. 2. Upper panel: The surface brightness profile of 3C 279 
(filled symbols) together with the model (solid line), the nucleus 
(dashed line) and the host galaxy (dot-dashed line). Lower panel: 
The surface brightness profiles of 3C 279 (solid line and filled 
symbols) and stars 1-4 in Fig.[T]relative to star 1. Each profile is 
traced out to a radius where the surface brightness can be deter- 
mined to better that 10% accuracy. Note that star 1 is represented 
by a horizontal line at SB - 0.0 in the lower panel. 

standard deviations above were chosen to mimic the range of 
PSF shapes seen in our i-band images and to roughly reproduce 
the residuals seen in the model subtracted image of 3C 279. 

We also included in the simulations the possibility that the 
host galaxy does not exactly follow the de Vaucouleurs profile 
with the slope of the surface brightness profile yS equal to 0.25, 
as often is observed (e.g. Nilsson et al.lll999l) . Therefore, when 
creating the host galaxy model the profile slope fi of the host 
galaxy was drawn from a Gaussian distribution with mean = 0.25 
and cr = 0.06. When performing the fits to the simulated images 
the j6 parameter was held constant at 0.25, however. 

After completing the simulations we computed the standard 
deviations of /core, /host and r^a to be 0.02 mag, 0.18 mag and 1.1 
arcsec, respectively. Taking into account the error in the calibra- 
tion we end up with /host = 18.4 + 0.3 and rgff - 2.7 + 1.1 arcsec 
for the host galaxy. The in the simulated fits had an average 
of 1.11 and standard deviation cr = 0.08. Given that the scatter 
in the simulations is significantly larger than the difference 
mx^ between the de Vaucouleurs and disk models, we conclude 
that we cannot claim either of the two models to be significantly 
better. 



4. Discussion 

Since most host galaxy parameters and correlations have been 
discussed in the R-band we first convert the host galaxy magni- 
tude from I-band to R-band. For the conversion we used A/ - 
0.06 for the galactic extinction (Schl egel et alJll998[) . Ki + e(I) 
= 0.4 for the K-correction and evolutionary correction, based on 



10g(MM,/Mo) 
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ref. 


8.912 


FWHM (Hp) - Up, 
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9.099 
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References: (1) iGu et ai ! ) ([200lh. (2) ICao & Jiarig ([2O02h. (3) 
Woo & Urrv (2002), (4) Liang & Liu (2003h. (5) IWang et al.l ( l2004h , (6) 
Bian & Zhag ( .2004.) . (7) .Liu et al, C2006) 



'Maraston! (Il998l) and lMaraston & ThomasI (l2000 l) and assuming 
galaxy formation at z - 2 an d finally R - 1 = 0.7 for an elliptical 
galaxy dFukugita et al.|[T995b . Using these values we derive R = 
18.7 for the host galaxy of 3C 279, which translates into Mr = 
-23.8 in the adopted cosmology. This value is in the bright end 
of the values found for the host galaxies of radio-loud quasars 
(-23.2 ± 0.3, Dunlop et al. 2003, transferred to the cosmology 
used h ere) and BL Lacertae objects (-22.9 + 0.5. ISbarufatti et al] 
|2005l). The effective radius is (17 + 7) kpc, which is fairly large 
compared to typical effective radii (~ 10 kpc) found for radio- 
loud quasars and BL Lacs, although the margin of error is also 
large. 

Next we estimate the central black hole mass of 3C 279 using 
the relationship between bulge luminosity and central black hole 
mass derived in Bentz et al. (2009) using the same cosmology as 
here: log(Mw,/lO^M0) = -0.02 + 0.8 log(Lbuige/lO'"L0), where 
Lbuige is measured in the V-band This relationship was derived 
using 26 Seyfert 1 galaxies and quasars at z < 0.3. Although the 
redshift of 3C 279 is outside this range, the host galaxy luminos- 
ity of 3C 279 is within the range of V-band bulge luminosities 
probed by this study 9 < log(Lbuige/i'0) < 11-5. 

As blazar host galaxies are expected to be bulge-dominated 
systems we use the derived R-band magnitude directly for com- 
puting the bulge luminosity. We first transferred the R-band 
lum inosity of the host g alaxy to the V-band using V - R = 
0.6 dFukugita et alJll995l) and used the above relation to derive 
log(M,,/,/Mo) = 8.9 for 3C 279. The scatter in the Lbuige - Mbi, 
relation is ~ 0.4 dex. The uncertainty arising from the host 
galaxy magnitude error is smaller (< 0.2 dex), so our esti- 
mate of the black hole mass has a total error of ~ 0.5 dex. 
Like the host galaxy luminosity, the derived black hole mass 
is in the upper end of the black hole mass distribution of 
FS RQs, but still consistent with the total observed range (see 
e.g. lGhiseUini & TavecchioH2 008). 

We have looked in the literature for other black hole mass 
determinations for 3C 279. The results of this search are sum- 
marized in Table |2] It is beyond the scope of this paper to dis- 
cuss in detail the various methods used to derive the values in 
Table |2] and their relative merits. Most of them have been de- 
termined with the empirical FWHM (H/j) - Lgpt correlation (e.g. 
IWandel et al. 1999) with different assumptions of BLR geome- 
try. There is a broad range of values found by different authors, 
but most of them are in the range 8.5 < log(Mw,/Mo) < 9.0, 
broadly consistent with the value found here. The two last val- 
ues in Table |2] deviate most from the general consensus and at 
least in Bian & Zhao (2004) (ref. 6 in table |2]i the high es- 
timate can be traced to a much broader reported line width 
than the other authors. 
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5. Summary 

The results of this paper can summarized as follows: 

(1) We have obtained a deep i-band image of 3C 279 which 
has enabled us to reveal the host galaxy for the first time. 

(2) The luminosity of the host galaxy (Mr = -23.8) is in the 
bright end of luminosities of other radio-loud quasar host galax- 
ies. The derived effective radius (17 + 7 kpc) is quite large com- 
pared to other radio-loud quasars, but the eff'ective radius is not 
very well constrained. 

(3) Using the empirical correlation between bulge luminos- 
ity and central black hole mass Mhh we estimate logCMfc/j/Mg) = 
8.9 + 0.5, broadly consistent with values obtained by photoion- 
ization methods. 
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